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O , Abstract 

I The decay 7r''e"*"e~ is being searched for as a signal of direct AS = 1 CP 

^P' violation. We provide a thorough updating of the analysis of the three components of 

the decay: 1) Direct CP violation, 2) CP violation through the mass matrix and 3) CP- 
conserving (two-photon) contributions. First the chiral calculation of the Ks — > 7T^e'^e~ 

■ rate, due to Ecker, Pich and de Rafael (EPR), is updated to include recent results on the 
. nonleptonic amplitude. Then we systematically explore the uncertainties in this method. 

J> I These appear to be so large that they will obscure the direct CP violation unless it is 

^rr ' possible to measure the Ks TT^e+e- rate. The CP -conserving amplitude remains 

. somewhat uncertain, but present indications are such that there may be a sizable CP- 

OO I violating asymmetry in the e^,e~ energies from the interference of the CP-conserving 

' and CP-violating amplitudes and this may potentially be useful in determining whether 

Q\ • direct CP violation is present. 

^: 

ji^' I. Introduction 

■ One of the goals of the next generation of rare kaon decay experiments is to attempt to 
^ ' observe CP violation in the decay Kl vr^e'^'e". This reaction is special because we expect 

^ ■ that direct CP violation (as opposed to the "mass matrix" CP violation already observed 
in the parameter e) may be the dominant component of the amplitude. This is in contrast 
with Kl — >■ Tin, where the direct effect is at most a few parts in a thousand of e. Direct CP 
violation distinguishes the Standard Model from "superweak"-type models [^. Moreover, the 
magnitude of the direct CP violation for this reaction is a precise prediction of the Standard 
Model, with very little hadronic uncertainty. In this article, we will update the analysis of the 
reaction Kl 7r^e~^e~ , attempting to understand if we can be certain that an experimental 
measurement is in fact a signal of direct CP violation. 

One difficulty is that there are three possible components to the decay amplitude: 1) 
A CP-conserving process which proceeds through two-photon exchanges, 2) a mass matrix 
CP-violating effect proportional to the known parameter e, and 3) the direct AS* = 1 CP- 
violating effect which we would like to observe. The existence of the first of these indicates 
that simply observing the total decay rate is not sufficient to unambiguously indicate the 
existence of CP violation. We need to either observe a truly CP-odd decay asymmetry, or else 
be confident on the basis of a theoretical calculation that the CP-conserving effect is safely 
smaller than the experimental signal. Unfortunately, the predictions in the literature for each 
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of the components listed above exhibits a range of values, including some estimates where all 
three are similar in magnitude. However, the quality of the theoretical treatment can improve 
with time, effort and further experimental input. We will try to assess the present and future 
uncertainties in the theoretical analysis. 

There remain significant experimental difficulties before it is possible to mount a search 
sensitive to a branching ratio of a few times 10~^^. We will assume that such a sensitivity is 
reached. At the same time, it is reasonable to assume that we will have improved experimental 
information on the related rate — > 7r''77, and that theoretical methods have provided a 
consistent phenomenology of this reaction. 

The KM parameters will be somewhat more fully constrained in the future, but hadronic 
matrix element uncertainties will prevent a precise determination of the parameters relevant 
for CP violation, at least until B meson CP violation has been extensively explored. With 
these expectations as our framework, will we be able to prove that the future experimental 
observation indicates the presence of direct CP violation? 

Our analysis indicates that one will not be able to prove the existence of direct CP violation 
from the branching ratio for — ^ 7r°e"'"e~ unless the decay rate for the related decay Ks — »• 
7r°e"'"e~ is also observed experimentally. This is yet more difficult than measuring the 
decay, and poses a problem for the program of finding direct CP violation. It is possible but 
not certain, that the electron charge asymmetry can resolve this issue and, when combined 
with the rate, signal direct CP violation. 

II. Overview of the analysis 

There is an extensive analysis associated with each of the three components of the decay 
amplitude which were listed in the introduction. We will devote separate sections of the 
paper to each of the major issues. The purpose of the present section is to highlight the main 
issues which are to be discussed more fully later, and to indicate how they fit together in an 
overall description of the decay process. 

In a way, the direct CP component is the simplest. The uncertainties are only in the basic 
parameters of the Standard Model, i.e., the mass of the top quark and the KM parameters. 
The relevant hadronic matrix element is reliably known. Unfortunately the extraction of KM 
elements has significant uncertainties, so that only a range of possible values can be given. 
This range corresponds to 7r°e"*"e~ branching ratios of a few times 10~^^. We discuss 

this range in Sec. III. 

The contribution of mass-matrix CP violation is more uncertain. The rate due to this 
source is given by the parameter e times the rate for Ks vr^e^e", so that the issue is the 
prediction of the CP-conserving Ks partial rate. Here the primary tool is chiral perturbation 
theory, with the pioneering treatment given by Ecker, Pich and de Rafael (EPR) @, |^. In 
Sec. IV, we update their analysis, under essentially the same assumptions. The main new 
ingredient is the inclusion of the results of the one loop analysis of nonleptonic decays, which 
decreases the overall strength of the weak K —>■ tt transition. This yields a change in the 
weak counterterms and a decrease of the rate. However, more important is an assessment of 
the uncertainties of such a calculation, which we describe in Sec. V. Any such calculation 
has a range of uncertainties, most of which we are able to estimate based on past experience 
with chiral calculations. We systematically discuss these. Unfortunately we find that one 
issue in particular has a devastating sensitivity on this mode. In their analysis, EPR made an 
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assumption which hes outside of chiral perturbation theory, that a certain weak counterterm 
satisfies W2 = 4:Lq where Lg is a known coefiicient in the QCD effective chiral Lagrangian. 
This assumption has no rigor, and the decay rate is very sensitive to the deviation W2 — 4:Lg. 
For any reasonable value of direct CP violation, there is an equally reasonable value of W2 
which can reproduce the corresponding Kl TT^e^e~ decay rate. Given a measurement, we 
will then be intrinsically unable to decide if it is evidence of a nonzero value of direct CP 
violation or merely measures a value for W2- It is this which indicates a need to measure the 
rate Kg — > 7r°e+e~. 

The third component is the CP-conserving amplitude which proceeds through the two- 
photon intermediate state Kl — > 7r°77,77 e^e~, described in Sec. VI. Here we must 
first understand the process 7r°77. This has been calculated in chiral perturbation 

theory at one loop order and has been measured experimentally. While the shape agrees with 
the chiral calculation, the theoretical rate misses by a factor of three. This has prompted 
some reanalyses of the theory of — >■ 7r°77, which we will take account of. However, the 
field has not reached a full conclusion on this mode, and it is clear that in the future the 
experimental and phenomenological status of this reaction will undoubtedly improve. We 
study how possible resolutions of these analyses will influence the —>■ 7r°77 decay rate. 
Ultimately this component should be satisfactorily understood. 

The ultimate problem is then our inability to distinguish, in a measurement of the 
Kl — ^ TT^^e+e^ decay rate, the direct CP violation from the mass matrix effect. It is pos- 
sible that the electron energy asymmetry may allow us to make this separation. The electron 
asymmetry comes from the interference of the CP-conserving two-photon process (even un- 
der the interchange of e+e~) and the CP-violating one-photon process (odd under the e'^e~ 
interchange). For many values of the presently favored parameter range, this asymmetry is 
very large i.e., of order 50%. In this case its measurement is not far more difficult than a good 
measurement of the rate. If we in fact are able to reach an understanding of the two-photon 
process, through future phenomenology and experiments on Kl — > 7r°77, then the asymmetry 
depends most critically on the CP-violating amphtude. If there is no direct CP violation, 
there is then a correlation between the Kl — > 7r°e"'"e~ decay rate and the electron asymmetry, 
parameterized by the unknown coefficient W2. As we detail in Sec. Vll, the presence of di- 
rect CP violation would upset this correlation, and in many cases would lead to a drastically 
different relative size of the asymmetry vs. decay rate, often even changing the sign of the 
asymmetry. Thus the asymmetry may be used to signal direct CP violation. Unfortunately 
this method is not foolproof. There exists combinations of values of W2 and KM angles for 
whom the distinction between direct and mass matrix CP violation are not so great and will 
be muddied by the inherent uncertainties in the theory. In Sec. VIII we explore the use of 
Kl — Ks interference to sort out the direct CP-violating amplitude. 

Overall, our reanalysis indicates that the demands on the experimental exploration of 
this reaction are quite severe. The simple observation of a few events will not be sufficient to 
indicate direct CP violation. The measurement of an electron asymmetry requires more events, 
and may or may not resolve the issue. Only the simultaneous measurement of Ks 7r°e"'"e~ 
allows a convincing proof of the existence of direct CP violation. 
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III. Direct CP violation 



Direct AS = 1 CP violation is manifested in the "penguin" reactions pictured in Fig. 1. 
The QCD short distance corrections to this mode have been thoroughly analyzed to next-to- 
leading order by Buras et al. 0, and we will use their results. The primary weak operator 
responsible for the transition have the form 

T^W = [C7vifJ')Q7V + C-jaQia] , (1) 

where 



Qjv = (sd)v~A{ee)v, 
QjA = {sd)v-A{ee)A- 



(2) 



The dominant contribution to the imaginary part of the coefficient Cn comes from the top 
quark, so that this is truly a short distance process. The coefficients have a CP-violating 
component 



Im Cji 

with the results of Ref. M yielding 



0.708 mt = 150, 
— = I 0.743 mt = 175, 
" ' 0.775 mt = 200, 



-Im f ^ ] y 



(3) 



yjA 

a 



-0.579 mt 
—0.736 mt 
-0.905 mt 



150, 
175, 
200, 



(4) 



with very little dependence on (the above is for Ajjg 

dependence on the low energy scale fi. 



0.3 GeV) and a negligible 




Figure 1: "Penguin" diagram inducing AS = 1 CP violation. 

The matrix element involved is well known via isospin symmetry from the charged current 
kaon decay, i.e.. 



(Ap') \dj,s\K%p)) = l^ip + p), + l^{p - p')„ (5) 

with 



A = (0.65 ± 0.005) W. (6) 

The form factor /_ does not contribute significantly to the decay because its effect is propor- 
tional to me- The decay rate is 



IniA, = lmV,,V* = \VJ\VJ sin 6 = A'X\ (7) 

where VJ^ = |V^^|sin5, and A,X,ri referred to the Wolfenstein parameterization of the KM 
matrix 0. This results in 

2.1 X 10-12 \ /T \ \ 2 f 150 GeV 

BR{Kl ^ 7r°e+e-)dir = { 2.75 x lO-^^ I x for = I 175 GeV . (8) 

3.6 X 10-12 J \10 / [200 GeV 

The dependence on the top quark mass will of course be removed by a convincing precise 
measurement of rrit. 

The KM parameter Vcb has the most favored values (including the recent CLEO data) 

0.036 ± 0.002 ± 0.002 HQET, 

_ ■ 0.036 ± 0.002 ± 0.003 B D*£iy models, , . 

~ ^ 0.039 ± 0.001 ± 0.004 B Xiu, 

0.036 ±0.003 (averaged). 

The element V^^ is measured by the inclusive decay B X^eu in the electron endpoint region. 
The two inclusive calculations available yield 

Kb _ f 0.082 ± 0.006 [ACCMM] [7], , . 

Kj, ~ j 0.074 ± 0.007 [RDB] §. ^ ' 

Models that calculate a set of exclusive decays {B —>■ Mev) can only be used to provide an 
upper bound on Kb since there are many final states (such as i? ^ vrTrez/ with tttt nonresonant) 
that are not calculated. These limits are 

r 0.12 [ISGW] [9], 
< <^ 0.087 [BSW] fig, (11) 

^•^b y 0.067 [KS] ini. 

We will use the former measurements to estimate 



K. 



0.078 ±0.007 ±0.010, (12) 
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with the first uncertainty experimental and the second theoreticaL In the Wolfenstein param- 
eterization of the KM matrix, the values of and imply 



A = 0.74 ±0.06, 
7]^ = 0.355 ±0.056. 



(13) 



Without any further analysis, these measurements imply an upper bound on ImAt: 



ImA* 



(1.0 ±0.3) X 10-Sin5. 



(14) 



A lower bound on this parameter can be found by consideration of the analysis of e. In the 
Wolfenstein parameterizations one has the approximate form 



|e| = (3.4 X W-^)A^r]BK 



l + 1.3A\l-p) 



rrit \ 



1.6' 



mw 



) 



(15) 



where Bk parameterizes the hadronic matrix element and is estimated to be in the range 
0.33 <Bk<1. Using (1 - p) < 1.4 and Bk < I one finds 



' 0.104 mt = 200 GeV, 
A'^T] > I 0.13 mt = 175 GeV, 
[ 0.16 mt = 150 GeV, 



(16) 



so that 



ImAt = A^X^rj > { 
This brackets the range 



' 5.3 X 10-5 mt = 200 GeV, 
6.8 X 10^5 mt = 175 GeV, 
8.2 X 10-5 mt = 150 GeV. 



(17) 



0.53 X 10^' < ImAt < 1-3 x 10"^ (18) 

Note that ImAt is positive. These constraints yield a decay rate from direct CP violation of 
magnitude 



1.01 ' 

1.25 i X 10-^^ < BR{Kl 
1.4 

Alternately, the "best" values 



7r°e+e )dir < < 



6.1 
4.6 
3.5 



> X 10 for mt = < 



f 200 
175 
150 



(19) 



ImAt 
mt 



1.0 X 10-^ 
175 GeV, 



(20) 



which we will take as our standard reference values, lead to a rate 
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BR{Kl 7r°e+e~)dir = 2.32 x lO"^^. (21) 
A more detailed analysis including a correlation between p and t] inherent in Eq. (15), as well 
as the use of B^B^^ mixing (which constrains (p — 1)^ + r]"^ as a function of /b) narrows 
the range only slightly because hadronic uncertainties dominate. 

IV. Revising the EPR analysis 

In this section, we review the formalism for analyzing mass matrix CP violation, first set forth 
by Ecker, Pich, de Rafael (EPR) [|, ||. This amounts to the prediction of the decay rate for 
Ks 7r°e+e~, since the mass matrix effect is defined by 

A{KL^7T'e+e-)\MM = eA{Ks ^ Tr'^e+e-), 

e = (2.258 X 10"^)e*^/^ (22) 

We then redo the results taking into account recent work on the nonleptonic kaon decays to 
one loop order. While this produces a significant numerical change, it is more important as a 
prelude to our subsequent analysis of uncertainties in the analysis. 

The prediction of Ks TT^e^e~ comes from a comparison with vr+e+e^, which 

contains many of the same ingredients. The reactions are displayed schematically in Figs. 2, 
3. In these diagrams the round circles represent the electromagnetic coupling while the square 
boxes indicate the action of the weak interaction. We know the electromagnetic interactions 
of pions and kaons from direct measurement. The weak K ^ it transition of Fig. 2a,b is 
known within some theoretical uncertainty from the use of chiral symmetry to relate it to 
K 2tt and if — > 37r. However, the weak Kn'y vertex is not known a priori and needs to be 
extracted from the analysis of 7T~^e~^e~ . 

The nonleptonic weak interactions are described by effective chiral Lagrangians, organized 
in an expansion in powers of the energy, or equivalently in numbers of derivatives and masses. 
At lowest order, called order E"^, the physical transitions are described by a unique Lagrangian 



C = GsTriXeD^UD^V^), 

U = exp ii ^ j, ^=1...8, (23) 

where 0"^ are the octet of pseudoscalar mesons (vr, K, rf). At next order, order i?^, the number 
of possible forms of Lagrangians is quite large, and has been categorized by Kambor, Missimer 



and Wyler [|1^. Not all of these contribute to K — 27r and K ^ -k^ but certain linear 
combinations do influence these amplitudes. The formalism of chiral perturbation theory 
dictate that when an analysis is carried out to order E"^, that one use Eq. (23) at tree level, 
in which case one obtains from the K it decay rate 

Gf 

Gs = 5 
g'r = 5.1. (24) 
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(c) 



Figure 2: Diagrams contributing to the reaction vr+e+e . Round circles represent the 

electromagnetic coupling while the square boxes indicate the action of the weak interaction. 




Figure 3: Same as in Fig. 2 for Kl 



Note that we neglect the CP violation in the nonleptonic amplitude, contained in Gg, 
because this is bounded to be tiny by the smallness of e'/e. 

In contrast, when evaluated at order E'^, one must include one loop diagrams in addition to 
the possible order E'^ Lagrangian. The loop diagrams involving tttt rescattering in the / = 
channel, K — > (7r7r)/=o — ^ (7r7r)/=o as pictured in Fig. 4, are quite large and are the major 
part of the order E^ analysis. While there is some ambiguity in the extraction of (see 
below), the enhancement from tttt rescattering leads to a smaller value of g^, with a good 



estimate being [jT3| 



gT = 4.3. (25) 

The K ^ IT amplitude used in Figs. 2, 3 does not have the enhancement from tttt rescattering, 
and is given in terms of g^ by 
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(26) 



We will explore further uncertainties in the K ^ tt vertex in the next section. 

To complete the diagrams of Figs. 2, 3 requires the electromagnetic vertices of kaons and 
pions. In chiral perturbation theory to order these are given by 



(7r+|J'^|7r+) 
{K+\J''\K+) 



1 + r 



TT 
TT 




967r2F2 



The first of these is known more fully from experiment, and has the form 



(tt+IJ'^Itt 



m = 730 MeV. 



Taylor expanding the latter form one determines Lq{ij, 
experimental charge radii 



(28) 

(7.4 ± 0.7) X 10-3. The 



(0.44 ±0.02) fm^ 
(0.34 ±0.05) fm^ 



(29) 



are compatible with this value. The final ingredient required for Figs. 2c, 3 is the weak 
photonic coupling. This includes both short distance and long distance physics, as illustrated 
in Fig. 5. While the short distance components have a reliable hadronic matrix element 
(it is due to the real parts of the coefficients discussed in the previous section), the QCD 
coefficient depends strongly on the low energy cutoff /i. In the full matrix element, this 
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dependence is canceled by a corresponding dependence on /x of the long distance physics. 
However, since the long distance physics is not calculable, we must attempt to determine this 
coupling phenomenologically. The innovation of EPR was to elucidate the possible forms that 
this coupling could take. They found that there were two possible chiral Lagrangians which 
could contribute to this process: 



w 



(30) 



In the presence of the short distance electroweak penguin effect due to Z exchange with an 
axial electron coupling, we need a third effective Lagrangian, not present in EPR, 



r' 



i27ra 



w 



The correspondence with the notation of the last section is 

3 1 Vya 



Im Wr, 



An \V ,V* I 

I ud usl 



a 



(31) 



(32) 



K 




TT 



Figure 5: Diagrams contributing to the short distance weak photonic coupling for K ire^e . 

Note that the labeling of takes into account the Lagrangians labeled by w^, W4 defined 
by EPR P, 1^, |l^ which contributed to other radiative K decays. In fact this form of is 
closely related to wi term since by using identities of the U matrix it can be shown that 



= {d,FnTr{QXe^M 
= {-e)e'^^eTr{QXQ^i'jCy) 

= e7^eTr(A6_i70, (33) 



10 



where in the second hne we have integrated by parts thus subsequently used the equation of 
motion so that 

^-^^F^^WiTriQXe-irCuJ^u) = ^^^'^^ wie-fi^eTr{X6-i7Cu) ■ (34) 
This aUows us to identify the short distance CP-violating part of wi. 

Im = ^ImAt. (35) 

The real parts of Wi, W2 contain long distance contributions and hence are not predictable 
by present techniques. They need to be extracted from experimental measurements. The 
other process available for this procedure is — > 7i~^e~^e~ (unless Ks 7r°e"'"e~ is measured 
in the future). Unfortunately one cannot fix both wi, W2 in this way, so that one must add a 
theoretical assumption in order to proceed. The relevant amplitudes are 

M(X+ ^ TT+e+e") = ^d+{p + p'Yu-f.v, 

4:71 

M{KL^n^e+e-) = -^{p + p'ru[dv7^. + dAW]v, (36) 

with 



d+ = W+ + (/)K{q^) + (t)n{q^), 

ds = Rews + 24>K{q^), 
167r2 



w+ = 



(w[ + 2wl-12Ll)--\og 
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w+ + IGtt^ {w; - 4LI) + - log 



m: 

1' 



f 

Jo 



dx 



and for the CP-violating Kl decay 



6 

1 7^x(l — x) 

mf 



m 



log 



2 ' 

K 

1 — — x) 
mf 



(37) 



dv 
dA 



^ds r — ^lm wi, 



ilm W5. 



(38) 



The goal of the search for direct CP violation is to separate the Im wi^s terms from the mass 
matrix effect eds- Note that in these expressions we have neglected the possible direct CP 
violation in the K ^ n transition (which is bounded to be very small by the measurement of 
e'/e) and the contribution of Re to CP-conserving decays (since Re W5 ~ Im W5 << Re Wi). 
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The EPR analysis of vr+e^e" uses the tree level value of g^, g^^"^ = 5.1. The decay 

rate is consistent with two values of Re w+, and a subsequent analysis of the decay spectrum 
favored the lower value for Re wi, i.e. Re W-^ = 1.16 ± 0.08 0. However, given that one 
is working to one loop order, it is more consistent to use the one loop value for g^, (73°°^ — 
4.3. Because of the presence of the Lg term, this is not just a rescaling of the value of 
An additional change that we make is to use the known full electromagnetic vertex in the 
pole diagrams, rather than just the first term in the expansion of the form factor. Note that 
because of the factor of A; ■ p in the weak matrix element, the only significant form factor is 
that of the pion in Fig. 2. This implies the replacement 

'Ff " (1 - qyml) ^^^^ 

in the formula for w^. The associated logarithm with fi ^ are much smaller than the 
Lg dependence and are not influenced much by this replacement. As a technical note, we 
comment that some potential modifications using a phenomenological pion form factor could 
lead to a lack of gauge invariance. By modifying the coefficient of a gauge invariant effective 
Lagrangian, we preserve the gauge invariant nature of the amplitude. With these changes, we 
find 



Re w;+ = 1.01 ±0.10. (40) 

[Without the second change, we would have had Re w+ = 1.33 ± 0.065]. This is illustrated in 
Fig. 6. 



BR(K"'^7T"'e'"e )xio'' 





Figure 6: The branching ratio BR{K^ vr+e^e^) is plotted against Re w+. The solid curves 
are obtained using the extreme values of the error intervals of \Vud\ = 0.9753 ± 0.0006, \Vus\ 
= 0.221 ± 0.003 and Lg{m^) = (7.4 ± 0.7) x lO"'^, while the dashed curve corresponds to the 
central values. The experimental value of the branching ratio ± its experimental errors are 
indicated by dashed horizontal lines. 

One cannot simply transfer this information to Ks or decays, because a different linear 
combination enters 
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1 Tfi^ 

ws = w+ + 167r\w2~AL,) + -log (41) 

mj^ 

However EPR deal with this problem by making the assumption that W2 = 4:Lg resulting in 
Re Wg = 0.73 ± 0.08 They note that this is an assumption which is not part of chiral 
perturbation theory, but do not explore the consequences if it is not correct. We will discuss 
this in the next section, finding a very strong sensitivity. At this stage we note that if one 
makes the assumption of W2 = 4:Lg, one obtains Re Wg = 0.58 ± 0.10 for our value of Re w+. 
At this value of Wg, the direct and mass matrix contributions are comparable: 

dv = eds zlm wi 

3 

^ e'''/^ (0.57 X 10"^) -il.O X 10"^ 
dA = il.O X 10-^ (42) 

when evaluated with rrit = 175 GeV, ImA^ = 10^^. This leads to a branching ratio 

BR{Kl 7r°e+e")MM = 0.37 x 10"^^ (43) 

if there is no direct CP violation {ImXt = 0) [EPR found BRmm = 1-5 x 10^^^ in this case 
[§],vs. 

BR{Kl 7r°e+e")cp = 1.78 x 10"^^ (44) 

for the full set of parameters given above in Eq. (37). The addition of mass matrix CP 
violation in this analysis led to a small decrease in the rate compared to the purely direct CP 
violation of the previous section, Eq. (21), because of the cancellation in the imaginary part 
of dv- However this may change if W2 ^ 4L9. 

One of the results that we will see in the next section is that the mass matrix contribution 
to the branching ratio is near a minimum value when W2 is close to ALg. For other values of 
W2, the rate can easily be an order of magnitude larger. Although our value for the branching 
ratio is a factor of four below that of EPR, both estimates are similar in saying that the mass 
matrix contribution will be small as long as W2 = 4:Lg. 

Using the tree level value for g^, gf'^^ = 5.1, the mass matrix contribution to the branching 
ratio is 

BR{Kl 7r°e+e")MM = 0.55 x 10"^^ (45) 



V. Uncertainties in mass matrix CP violation 

Our goal in this section is to assess how well we understand the prediction for Ks — * 7c^e~^e~ . 
The most important effect will be discussed in section (c) below, but we proceed systematically 
to discuss even contributions which have less uncertainty. 

a) Purely electromagnetic vertices. 
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The electromagnetic vertices enter in diagrams 2, 3a,b. The uncertainty here is in the 
choice of whether to use the chiral expansion of the form factor truncated at order g^, Eq. 
(27), or the full dependence of the monopole form factor, Eqs. (28), (39). The first 
choice is natural when one is working to a given order in the chiral energy expansion, 
but the latter choice clearly includes more of the physics which is known about the 
electromagnetic vertex. Note that it is essentially only the pion form factor which is 
relevant, because the diagram involving the kaon form factor. Fig. 2b, is suppressed by 
a factor of m^/m|- with respect to Fig. 2a because of the momentum dependence of the 
weak K ^ 71 transition. The use of the full form factor produces a modest variation in 
the value of w+ (i.e. Re w+ = 1.01 instead of Re = 1.33). Because of cancellations 
in the amplitude this provokes a more extreme variation on the decay rate. The 
lowest order chiral vertex, Eq. (27), produces a decay rate 

BR(Kl 7r°e+e-)MM = 1-96 x 10"^^ (46) 

instead of the result of Eq. (43). [We note that if we had also modified W2 in the same 
way as ALg as in Eq. (39) we would have a BR{Kl — > 7r'^e+e~)MM = 1-33 x lO""-*^^.] 
While we feel that it is good physics to use the full electromagnetic form factor, one 
could also interpret these results as an uncertainty in the analysis due to higher order 
terms in q^, with that uncertainty being of order 2 x 10~^^. 

b) The weak X — > tt vertex. 

Wc have already given one indication of the sensitivity of the result to the size of the 
K ^ n transition. Under otherwise identical assumptions, gl'^'^'^ = 5.1 yielded the rate 
in Eq. (45), while 513°°'' = 4.3 produced the result of Eq. (43). These modifications 
to also include corresponding modifications to the Ktt^ vertex required by gauge 
invariance. This is automatically maintained, however, by the use of gauge invariant 
effective Lagrangians. In order to appreciate that this change in is not the only 
uncertainties in the Kn amplitude, one needs to understand a bit more about the chiral 
phenomenology of X — > 27r and K 3n. 

Chiral symmetry relates processes with different numbers of pions, such as i^T — > tt vs. 
K — > Sir. The predictions are compactly contained in the chiral Lagrangians, but can 
also be obtained using the soft pion theorems, which was the methodology used in the 
1960's. The only advantage of the latter technique is that it relies only on chiral SU(2) 
while modern chiral Lagrangian analyses have always involved chiral SU(3) symmetry. 
[Presumably the latter could be reformulated in chiral SU(2), but no one has yet done 
this.] The soft pion analysis indicates that one obtains the same relation (up to terms 
of order m^) between X — > 27r and K ^ Sir for any Lagrangian which survives in any 
soft pion limit of — 37r (i.e., pj — » 0). The only Lagrangians which do not survive in 
any soft pion limit involve four separate derivatives on the four fields of K ^ Sir, e.g. 

>Cquartic = j^Tt (X^D D D'^U D'' U^) . (47) 

This Lagrangian yields a matrix element proportional to {k ■ Pi){p2 • Ps) which clearly 
vanishes as any — > 0. In contrast most of the order Lagrangians do not vanish in 
all soft pion limits. An example is 
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= ||Tr (u^XeD^DMD^'U^D''U) , (48) 



which then yields the same relations of K 3tt and K —>■ 27c as does the lowest order 
result C given in Eq. (23). [This phenomenon is explained in more detail in Ref. |T3[.] 



Since the only inputs to the chiral phenomenology are the amplitude for K ^ 27t and 
K — *■ Stt, it follows that one cannot distinguish a combination C + C from a Lagrangian 
involving C only. 

However, when we discuss the K —>■ tt vertex, there is a distinction between these various 
Lagrangians. For example C in Eq. (48) involves a minimum of three meson fields, and 
hence contributes to K —>■ 27T and K ^ in but not at all to ^ tt. In contrast, 
the lowest order Lagrangian Eq. (23) contributes to all of K —>■ tTjK ^ 27r, K Svr. 
Since known phenomenology cannot distinguish between linear combinations of C + C, 
this manifests itself in an uncertainty in the K ^ tt vertex. Since the higher order 
Lagrangians are known to make a 25% difference in relation between K ^ 2i[ and 
K — s> Stt, it would be unreasonable to take this uncertainty in i^T — > vr to be any less 
than 25-30%. Using = 4.3 ■ (1 ± 30%) yields a range 

BR{Kl 7r°e+e-)MM ~ 0(10"^^) ^ L5 x 10"^^ (49) 

This again indicates that the analysis has significant cancellations present and that 
modest variations in the analysis can lead to uncertainties of order 1.5x10^^^. 

c) The Kfi'y vertex. 

The parameter w+ was determined from the analysis of — > 7r+e+e^. This arrange- 
ment has some intrinsic uncertainty because it was performed to a given order in the 
chiral energy expansion. It would be reasonable to take this uncertainty at 30%. 

However, this uncertainty is dwarfed by the errors introduced by the assumption of 
W2 = 4L9. There is nothing about chiral symmetry which forces such a relation. For 
example, if the process of Fig. 7 were to contribute to the weak coupling, the relation 
W2 = 4L9 would not occur except for a special value of the K ^ ai, fi amplitude. This 
is highly unlikely, and we can easily accept W2 — 4:Lg 7^ 0. The crucial distinction here 
is between models and rigorous theory. Chiral perturbation theory is a rigorous method 
which expresses true relationships in QCD to a given order in the energy expansion. 
However an assumption such as W2 = 4:Lg may be true or false in a way that we cannot 
decide based on QCD. It may occur within some models, yet we have no guidance 
as to whether it is correct in nature. We cannot base something as important as the 
observation of direct CP violation on something as flimsy as a model. 

Unfortunately the decay rate for Ks vr^e'^e" depends very strongly on the value of 
W2. This strong dependence was observed by Littenberg and Valencia |T^. If we were 



to chose W2 = 0, the rate would be two orders of magnitude larger. In Fig. 8 we plot the 
branching ratio for the mass matrix contribution to — > vr^e^e" vs. W2- We see that 
reasonable values of W2, we get a wide range of values of the branching ratio. Conversely, 
a measured value of BR{Kl 7r°e"^e~) in the range of 10~^^ —>■ few x 10~^^ could be 
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Figure 7: Diagram contributing to the weak photon couphng. The notation for the vertices 
is as in Fig. 2. 

interpreted in terms of a reasonable value of W2- This is then an enormous uncertainty 
in the mass matrix contributions to Tr^e+e". 

This uncertainty could be removed if one measured the rate of Ks vr'^e+e". The mass 
matrix contribution would then be known, see Eq. (22). However this task is not easy 
experimentally. 




Figure 8: The mass matrix contribution for the branching ratio BR{Kl 7r^e~^e )mm is 
plotted against W2. The convention for the sohd and dashed curves is the same as for Fig. 6. 



VI. The CP-conserving amplitude 

The Kl Ti^e^e^ transition can also take place through a CP-conserving two-photon inter- 
mediate state. If we ignore the electron mass, the form of the amplitude will be 

M{Kl ^ 7r°e+e-)cpc = G^a^Kp ■ {K+ - K-){p + p'Tu-f^^v, (50) 

where X is a form factor and the extra antisymmetry under ke+ ^ kg- is a reflection of the 
properties under a CP transformation. In order to calculate this we need to understand the 
Kl — > 7r°77 transition first. 
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We are fortunate that Kl — >■ 7r'^77 is accessible to present experiments, and significant 
phenomenology has been performed on this reaction. We will utilize the work of Cohen, 
Ecker and Pich Jl^ as representative of present work, with the understanding that future 
experimental and theoretical work will clarify the analysis considerably. The Kl —>■ 77^77 
amplitude has two form factors A, B, defined via 



M (Kl 7r°77) 
B 



^e/,(gi)e^(g2) 



tJ.U\ 



+2^- {k ■ qiqi^p" + p ■ q2qiP^ - qi ■ q2P^p'' - g'^'p ■ qiP ■ 92 ) 



m 



K 



(51) 



When the photons couple to e"'"e~, as in Fig. 9, it is well known that the A amplitude 
contributes to — ^ 7r°e"'"e~ only proportional to mg, which is a small effect which we will 
drop. It is the B amplitude which is important for the e^e~ final state. The authors use a 
representation which fits the known Kl tt+tt^tt^ amplitude in a dispersive treatment of 
Kl 71^77 and find 



B{x) 

X 
C2 

F{x) 



[1 4 

C2 -F{x) + -(5 - 2x) 

I Ob O 



+ R{x] 



2 , ml] 
6 nip ) 



R{x) 




X > 1, 



X < 1, 



+ m 



X > 1. 



(52) 



The most important ingredient above is ay which is an unknown parameter representing the 
vector meson exchange contributions to the B amplitude. A fit to the decay rate (using 
► 7r%7) 



BR{Kl 



1.7 X lo- 



ll 



and 77 spectrum in K^ — > 7r°77 indicates a value 
around ay = —0.96. This parameter was very important in increasing the chiral prediction 
of the decay rate to be in agreement with experiment. We have explored the sensitivity of 
this parameter to changes in the analysis and have found that 25% changes in the dispersive 
treatment lead to a factor of 2 change in ay, so that this value is still quite uncertain. 



Several authors [|14| , |17[ |l^ have calculated the contribution of the on-shell two-photon 
intermediate state to K^ — ^ Tt'^e^e" . Although this is sometimes referred to as the absorption 
contribution, it is not the full absorption part since there is a further cut due to on-shell pions. 
Besides this, the full CP-conserving amplitude also receives contribution form the dispersive 
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Figure 9: CP-nonviolating diagram involving two photons coupling to e^e . The notation for 
the vertices is as in Fig. 2. 



part of the amplitude, with off-shell photons (and pions). The calculation of this is complicated 
by the sensitivity of the loop integral to high momentum, as the Feynman diagram of Fig. 9 
will diverge if we treat the B amplitude as a constant in and g|. However the remedy to 
this is well known; the couplings of virtual photons to hadrons is governed by form factors 
which lead to suppression of the couplings at high q^. We will include an estimate of these 
form factors and this will allow us to calculate the dispersive component of the CP-conserving 
amplitude. 

The two-photon loop integral in the limit me — >^ is given by 



M{Kl ^ 7r°e+e-)cpc 



Gga f dH B{k + k')F{e + k)F{£-k') 



(27r) 



e^{i + ky{e-k') 



X 



{^[i^ + p ■ {k - k') - p ■ ££ ■ {k - k') + 2p ■ ek ■ k' - p ■ ke ■ k' - p ■ k'i ■ k] + 

/ [{P -ly+p-lp-ik- k')] } . 



(53) 



Here F(g^) are the form factors for the virtual photon couplings. The structure above is 
certainly an approximation, as in general the virtual photon dependence need not be only an 
overall factor of F[q^). However the above form would be sufficient to capture the kinematic 
variation if only one photon is off-shell (given an appropriate F{q^)). Since we only need a 
minor form factor suppression to tame the logarithmic divergences, we feel that this structure 
will be sufficient for our estimate. We choose 



-m 



V 



(g2 - rriy) ' 



(54) 



which is a good representation of almost any mesonic form factor, with my ~ mp. Neglecting 
terms which are suppressed by powers of l/rrip, we find the amplitude of Eq. (50) above, with 



K 



Six) 



167r^m|- 




(55) 



where s — {k^- + k^+Y ■ The log factor is of course expected, since the photon "absorptive" 
part comes from the expansion log(— s) = logs -|- ivr. 

This representation of the amplitude leads to a CP-conserving branching ratio of 



18 



BR{Ki 



e-)cpc = 4.89 X 10 



-12 



(56) 



for av = —0.96. More generally we show the CP-conserving branching ratio vs. ay in Fig. 
10. Note that while for many values of ay the CP-conserving rate is small compared to 
the CP-violating rate of previous sections, these two rates are comparable for some range of 
parameters. Note that there is no interference in the rate between the CP-conserving and 
violating components, so that the rates just add, as shown in Fig. 11. However, there can be 
a CP-odd asymmetry in the electron positron energies, which we turn to in the next section. 



BR(K°^7r°e"'e")xiO 



BR(K°^7r°e*e-)xi0' 




Figure 10: The CP-conserving branching 
ratio BR{Kl — > 7r°e"^e")cpc is plotted 
against ay. The convention for the solid 
and dashed curves is the same as for Fig. 
6. The assumed value for ay is indicated 
by the dashed vertical line. 



Figure 11: The complete branching ratio 
BR{Kl — > 7r°e+e~) is plotted against W2- 
The convention for the solid and dashed 
curves is the same as for Fig. 6. ay — 
—0.96 is assumed. 



VII. The Electron Energy Asymmetry 

When both CP-violating and CP-conserving amplitudes contribute to the decay, there will be 
an asymmetry in the electron positron energy distribution 

^ ^ N{E^ > E_) - N{E^ < E_) 

N{E+ > E_) + N{E+ <E_)' ^ ' 

This will be quite large if the two amplitudes are comparable. In contrast to the overall 
decay rate, this asymmetry is an unambiguous signal of CP violation. It may be even more 
useful if it can be used to prove the existence of direct CP violation. This can occur because 
the asymmetry is sensitive to the phase of the CP-violating amplitude, and mass matrix CP 
violation has a unique phase (that of e), while direct CP violation will in general have a 
different phase. 

For the electron energy asymmetry to be useful as a diagnostic of the form of CP viola- 
tion, the CP-conserving two-photon amplitudes must be well known. As we discussed in the 
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previous section, it is reasonable to expect that this will be true in the future after further 
phenomenology of n'^'j'j. For illustrative purposes, we will use ay = —0.96. The analy- 

sis below will need to be redone in the future if this value of ay changes significantly. However 
the pattern of the analysis and the general conclusions will be valid for a wide range of values 
of av- 

The amplitudes involved have been given in the previous section. Note that the two-photon 
amplitude has both real and imaginary parts, the mass matrix CP-violating amplitude has a 
phase of 45°, while the direct component is purely imaginary. The asymmetry is proportional 
to the imaginary part of B times the real part of dy minus the product of the real part of B 
and imaginary part of dy- The asymmetry may be defined in differential form: 

— (E+>E_) = I dy- 



dz Jo dydz 

— {E+<E_) = I , ^dy 



dz J-iA^(i,z,r2) dydz^ 

A5(l,z,r2) = 1 + z'^ + -2z -2r'^ -2r 



r 



A{z) 



rriK 

dT{E+ > E_)/dz - dT{E+ < EJ)/dz 
dT{E+ > E_)/dz + dV{E+ < E_)/dz' 



In Fig. 12 we plot the differential asymmetry for several values of W2 in the case when there 
is no direct CP violation. Fig. 13 gives the same information for ImAf = 10~^. We see that 
the asymmetry is sizable for many values of W2 and that it depends significantly on direct 
CP violation, even possibly changing sign. The integrated asymmetry is plotted versus W2 in 
Figs. 14, 15. 

Both the decay rate and the asymmetry depend on Wg and this forms the main uncer- 
tainty in the analysis. However, it is possible to remove this uncertainty by measuring both 
observables. For a given value of the CP-conserving amplitude, there is a strict correlation 
between the two. In Fig. 16 we plot the values of the branching ratio and the integrated 
asymmetry A as one varies W2- We see that the curves with and without direct CP violation 
are well separated for much of the range. To the extent that we understand the CP-conserving 
amplitude, this can be used as a diagnostic test for direct CP violation. 



VIII. Time dependent interference of Ki^ Tr^e^e 

Littenberg |^ has suggested a time dependent analysis of a state starting out as a in 
order to extract maximal information about the decay mechanism. This is far more demanding 
experimentally then simply measuring Kl or Ks decays separately. However we analyze this 
technique in order to assess its usefulness. 

A state that at t = is a K*^ will evolve into a mixture of Kl and Ks'- 

\K\t)) = j^^^ {e-^^^* [\K,) + e\K2)] + e"^^-* [jK^) + e\K,)]] , (59) 
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Figure 12: The differential asymmetry Figure 13: Same as Fig. 12 for ImAj — 

A{Kl n^e'^e') is plotted against z = 10~^. 

{k + k'y/mji for ^2 = 1 xlO~^ (curve 1), 

2 xlO"^ (curve 2), . . ., 5 xlO"^ (curve 

5), in tlie case wlien tliere is no direct CP 

violation, ay = —0.96 is assumed. 




Figure 14: The integrated asymmetry is 
plotted vs. W2 in the case when there is no 
direct CP violation. The convention for 
the sohd and dashed curves is the same 
as for Fig. 6. ay = —0.96 is assumed. 
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Figure 15: Same as Fig. 14 for ImAt = 

10-1 
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where 

Hj=m,-t^, J = S,L. (60) 

Ignoring small effects such as second order CP violation and direct CP-violating components 
in the parameter e, we then have a time development proportional to 



|(7r°e+e-|7i|i^°(i))p 



^ { \Asf e-^'' + \eAs + A^ir + ^cpc|' e"^^* 



+ 2Re 



{eAs + Adir + Acpc) A* e-^("^--"^«)*l e"^^^ \ . (61) 



Here A^ir, eAs and Acpc are the amplitudes analyzed in Sees. Ill, IV and VI respectively. 
Measurements at early time {t ^ l/P^) determine T{Ks — > vr'^e+e^) while at the late time 
(t ^ l/r^) one observes T{Kl 7r^e~^e~). These contain the information described in 
preceding sections. However in the interference region t = O {l/{mL — ms)) ~ 0(l/r_g), we 
obtain extra information about the separate contributions to the decay amplitude. 

In Fig. 17 we show the time dependent signal for the cases of pure mass matrix CP 
violation and the addition of direct CP violation, using the analysis of the previous sections 
with W2 = 4:Lg. We see that the shape of the interference region does differentiate these two 
cases, but that for the case studied the dependence on direct CP violation is not so large as 
to allow an easy experimental determination. 
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Figure 16: The integrated asymmetry 
A is plotted vs. the branching ratio 
BR{Kl — > 7r°e"'"e~) as one varies W2- 
ay = —0.96 is assumed. The solid curve 
describes the case when there is no direct 
CP violation, the dashed curve is for ImAt 
= 10-^ 



Figure 17: The normalized time distribu- 
tion of 7r^e^e~ decays is plotted 
vs. time t in usees, ay = —0.96 is as- 
sumed. The solid curve describes the case 
when there is no direct CP violation, the 
dashed curve is for ImAt = 10~^. 
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IX. Summary 



Because there are three possible contributions to Kl — > vr^e+e^, the analysis has multiple 
issues which theory must address. We have provided an updated analysis of all of the compo- 
nents. The goal of identifying direct CP violation will not be easily accomplished. The decay 
rate by itself suffers from a severe uncertainty in the analysis of the mass matrix contribu- 
tions. In the chiral analysis there is a free parameter, ^2, which is not fixed experimentally, and 
which has a strong influence on the decay rate. Measuring the related rate of Ks ~^ n^e~^e~ 
would determine this parameter and will likely allow us to determine whether or not direct 
CP violation is present. 

Alternatively if the Ks branching ratio is not measured it may be possible to signal direct 
CP violation using the asymmetry in the e"^, e" energies. The direct and mass matrix CP 
violations have different phases, and the asymmetry is sensitive to their difference. There 
exist some combinations of the parameters where there remains some ambiguity, but for 
sizeable portions of the parameter space direct CP violation can be signaled by a simultaneous 
measurement of the decay rate and the energy asymmetry, as illustrated in Fig. 16. 

Progress in the theoretical analysis is possible, especially in the CP-conserving amplitude 
which proceeds through the two-photon intermediate state. Here both a better phenomeno- 
logical understanding of the related decay — ^ 7r°77, and a better theoretical treatment of 
the dispersive contribution should be possible in the near future. 
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